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PROBABILITY DISTRIBUTIONS FOR THE ERRORS IN THE 
PARAMETERS OF NEAR-EARTH CIRCULAR ORBITS 
WITH APOLLO APPLICATIONS 


C. W. Murray, Jr. 


ABSTRACT 

A statistical technique is described for analyzing the effects of 
insertion errors (actual minus nominal) on the parameters of near- 
earth circular orbits. Probability distributions are obtained 
numerically using Keplerian two-body equations of motion, and it is 
shown that the errors in perigee, apogee, and eccentricity are non- 
Gaussian with non- zero means. The technique is applied to a nom- 
inal Apollo earth parking orbit. The analysis indicates that the 
actual near-earth parking orbit will be close to the nominal 100 
n.mi. circular orbit. For example, there is a 99.5% probability that 
the actual perigee height will exceed 97.5 n.mi., and a 99 .5% prob- 
ability that the actual apogee height will be less than 102.4 n.mi. 
Insertion errors (actual minus nominal) and insertion ship tracking 
errors (calculated minus actual) are combined as a further example. 
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PARAMETERS OF NEAR-EARTH CIRCULAR ORBITS 
WITH APOLLO APPLICATIONS 

C. W. Murray, Jr. 


SUMMARY 

A statistical technique is described for analyzing the effects of insertion 
errors (actual minus nominal) on the parameters of near-earth circular orbits. 
The error in the state vector (position and velocity) at insertion is assumed to 
have a multivariate normal distribution with a zero mean and given covariance 
matrix. Probability distributions for some of the errors in the parameters are 
obtained numerically using Keplerian two-body equations of motion, and it is 
shown that the error in perigee, apogee, eccentricity, and the angle between the 
actual position vector and the nominal position vector at insertion are non- 
Gaussian with non- zero means. 

The technique is applied to a nominal Apollo earth parking orbit using an 
expected covariance matrix of insertion errors typical of the performance of the 
Saturn V Launch Vehicle at insertion into a near-earth nominally circular park- 
ing orbit of 100 n.mi. The analysis indicates that the actual parking orbit will be 
close to the nominal 100 n.mi. orbit. For example, there is a 99.5% probability 
that the actual perigee height will exceed 97.5 n.mi., and also a 99.5% probability 
that the actual apogee height will be less than 102.4 n.mi. 

Insertion errors (actual minus nominal) and insertion ship tracking errors 
(measured or calculated minus actual) are combined as a further example. Re- 
sults indicate that under worst case conditions for the coefficients of correlation 
between the tracking errors in insertion height, speed, and flight path angle (co- 
efficients of +0.9), there is a 90% probability that the calculated perigee height 
will exceed 91 n.mi. The three sigma values for the insertion ship tracking 
errors in insertion height, speed, and flight path angle used in the analysis are 
(Reference 1): 


3°Ar = 2.4 n.mi. (4.44 km) 

3°A v -*-6 ft/sec (4.87 m/sec) 
3 < 7 A r , 0.16° (2.79 mrad). 


iv 


CONTENTS 


Page 

ABSTRACT iii 

SUMMARY iv 

INTRODUCTION 1 

1. General 2 

2. The(U°V°W°) Coordinate System 3 

3. Expressions for Some of the Errors in the Parameters 5 

4. The Variances and Covariances Between the Errors in Insertion 

Radius, Speed and Flight Path Angle 11 

5. The Cumulative Distribution Functions of the Errors in the 

Parameters 12 

6. An Example With Application to a Nominal Apollo Earth Parking 

Orbit 16 

CONCLUSIONS 29 

ACKNOWLEDGEMENTS 30 

References 30 


v 


ILLUSTRATIONS 


Figure Page 


1 The (U° V° W°) Coordinate System 4 

2 The Cumulative Distribution Function of the Error in 

Insertion Radius AR 0 19 

3 The Cumulative Distribution Function of the Error in 

Insertion Speed AV 0 20 

4 The Cumulative Distribution Function of the Error in 

Insertion Flight Path Angle AT 0 20 

5 The Cumulative Distribution Function of the Error in 

Energy Per Unit Mass AC 3 21 

6 The Cumulative Distribution Function of the Error in 

Semi- major Axis AA 21 

7 The Cumulative Distribution Function of the Angle Between 
the Actual Position Vector and the Nominal Position 

Vector at Insertion B 22 

8 The Cumulative Distribution Function of the Error in 

Eccentricity AE 22 

9 The Cumulative Distribution Function of the Error in 

Perigee Radius AR p 23 

10 The Cumulative Distribution Function of the Error in 

Apogee Radius AR a 23 

11 Comparison of the Cumulative Distribution Function of the 
Error in Perigee With a Cumulative Normal Distribution 

Having the Same Mean and Standard Deviation 25 

12 Probability That Calculated Perigee Height Exceeds r ... 28 

P t 



LIST OF TABLES 


Table Page 

1 Expected Saturn V Insertion Covariance Matrix 16 

2 Correlation Matrix for Matrix in Table 1 17 

3 Covariance Matrix £ 0 of AR 0 , AV 0 . and AT 0 18 

4 Normalized Covariance Matrix of AR Q , AV 0 , and AT 0 19 

5 Lower and Upper 99.5% Probability Points for the Errors 

in the Parameters of the Orbit 24 

6 Covariance Matrix of AR X , AV t , and AT X 27 

7 Covariance Matrix 2 of (AR q +ARJ), (AV 0 +AVJ), and 

(ATo+AT,) 28 



PROBABILITY DISTRIBUTIONS FOR THE ERRORS IN THE 
PARAMETERS OF NEAR-EARTH CIRCULAR ORBITS 
WITH APOLLO APPLICATIONS 


INTRODUCTION 

The near-earth orbit of a space vehicle can be determined from the position 
and velocity vectors of the vehicle at insertion using Keplerian two- body equa- 
tions of motion. Due to insertion errors (errors in the xyz components of these 
vectors within an inertial Cartesian coordinate system at the time of insertion), 
the actual orbit will deviate from the desired nominal. 

The purpose of this report is two-fold: (1) to study the probability distribu- 
tions of some of the errors in the parameters of near-earth circular orbits 
assuming the insertion errors are correlated Gaussian errors with zero means; 
(2) to demonstrate a statistical technique using numerical integration for deter- 
mining the probability distributions of non-Gaussian errors. 

By stating that the errors in the Cartesian coordinates of the vehicle's 
position and velocity vectors at insertion are Gaussian, we are making a more 
basic assumption— that these particular errors can be represented as linear 
combinations of a number of error sources within the guidance and control sys- 
tem of the launch vehicle (gyro drift, accelerometer errors, etc.) which are 
Gaussian (but not necessarily independent).* The assumption of linearity is 
reasonable if second order effects (and higher) of these error sources are quite 


*A numerical procedure for determining the coefficients of the error sources in the linear expression 
for each of the errors in the state (Ax, Ay, Az, Ax, Ay, Az) is to take each error source one at a 
time, and using its one sigma value, determine the error in the six-dimensional state vector 
through a digital computer powered flight program which simulates the guidance system in the 
launch vehicle. If there are N sources of error, N + 1 computer simulations would have to be 
made (one run for the nominal state). In matrix notation, if AX is a (6 x 1) vector representing 
the error in the state (position and velocity) of the vehicle at insertion, Q is a (6 x N) matrix, 
and AW is an (N x 1) vector of error sources, 

AX = QAW 

and the above procedure is equivalent to determining one column of the Q matrix at each step. 
Then, the covariance matrix of the error in the state will be given by: 

M = QPQ t 

where P is the correlation matrix (l’s down the main diagonal and the coefficients of correlation 
as off-diagonal elements) of the (N x 1) error vector. If AW has mean vector £ (error sources 
with biases), then AX will have mean vector Usually these error sources are taken to be 
independent and to have zero mean. 
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small compared to first order effects (certainly this would be true in a well- 
designed guidance system). The Gaussian assumption is reasonable if one con- 
siders these error sources to be unpredictable— constant for one flight but varying 
in a random fashion over many flights in accordance with a Gaussian distribution. 
Even if this is not exactly true, the assumption that the error in the state (posi- 
tion and velocity) vector at insertion has a multivariate normal distribution with 
an associated covariance matrix and zero mean provides us with a working tool 
to comparatively analyze the performance of one space vehicle with another and 
to obtain bounds on the errors in the parameters of the orbit for a particular 
mission. 

It will be seen that some of the errors in the parameters of the orbit can be 
expressed to a good approximation as linear combinations of the insertion errors, 
and therefore, will be normally distributed (Gaussian). Others cannot, and the 
probability distributions of these will be non-Gaussian, for which a numerical 
technique is necessary to obtain the distribution. 

An example will be given which applies the technique to a nominal Apollo 
earth parking orbit. 

Determination of the probability distribution of a random variable is im- 
portant for two reasons: (1) It gives more information than just the mean and 
standard deviation, in particular, for non-Gaussian distributions; (2) It may be 
important for orbital parameters having critical upper or lower bounds such as 
perigee height in the case of the Apollo parking orbit. (Reference 2). 


1. GENERAL 

Let X represent the nominal state vector (position and velocity) of a space 
vehicle (within an inertial coordinate system) at insertion into a near-earth 
orbit, and let 77 be one of the parameters of the orbit. Then, using Keplerian 
two-body equations of motion, we can express 77 as a function of X 


V ~ t?(X) 


( 1 ) 


Due to insertion errors, the actual state X p will differ from the nominal 
state X by an amount Ax 


X p - X + Ax 


( 2 ) 


2 



We may therefore write an error At? in the parameter 7] as 


At ? - ^( X p) " _ t](X + Ax) - 7](X) (3) 

Let AX be a random vector representing the error in the state vector at 
insertion. 

Then 


AH = ?](X + AX) - t]( X) 


(4) 


is a random variable* representing the error in rj. 

In this analysis we will assume that the nominal orbit is circular and that 
AX has a multivariate normal distribution with zero mean and a given covariance 
matrix. 

It will be seen that some of the errors in the parameters of the orbit can be 
expressed to a very good approximation as linear combinations of the components 
of AX, and, therefore, will be normally distributed. Other errors (e.g., perigee 
error and apogee error) cannot be expressed as linear combinations, and for 
these a numerical technique must be used to determine the distribution. 


2. THE (U°V°W 0 ) COORDINATE SYSTEM 

The ( U° V° W°) system is a right-handed inertial Cartesian Coordinate Sys- 
tem frequently used to define the position and velocity of a space vehicle at 
insertion. It may be seen in Figure 1. 

The U° vector is a unit vector in the R (position vector) direction. The W° 
vector is a unit vector in the (R x V) direction (V is the velocity vector), and the 
V° vector is a unit vector in the (w° x U°) direction. 


*lt should be noted that Ax is a value which the random vector AX can take. Also, At] is a value 
which the random variable AH can take. 
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Figure 1. The (U° VO WO) Coordinate System 


In the (U° V° W°) System the R and V vectors for a nominally circular orbit 


are 


R = 


ro 

0 

0 _ 



( 5 ) 


( 6 ) 


c 
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where r Q is the radius at insertion (radius of the orbit) and /u. is the gravitational 
constant. The perturbed (or actual) position and velocity vectors are then 


R„ 


r o + Ax i 
Ax 2 

Ax, 


Ax. 

4 

V 0 + Ax 5 

Ax. 


(?) 


where v 0 - y'/x/r 0 . 


3. EXPRESSIONS FOR SOME OF THE ERRORS IN THE PARAMETERS* 


Insertion Radius Error 


The error in insertion radius is 


Ar o = r 0 (X+Ax) - r Q = /(r fl +AxJ 2 + (Ax 2 ) 2 + (Ax 3 )' 


r 


= ( r o +Ax i) 


1 + 


i\ K ) 2 + ( a * 3 ) : 




Ax l) : 


" r- r. 


Ax, 


neglecting 2 nd order terms and higher. Or, in random variable notation 


AR q % AXj 


( 8 ) 


* Thro ug hout this Section capital letters will denote random variables. Bold letters represent 
vectors and matrices. 
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Insertion Speed Error 


The error in insertion speed is 


Av o = v o (* + Aj O - v 0 = ]/(Ax 4 ) 2 + (v 0 +Ax 5 ) 2 + ( Ax 6 ) 2 v o 


r 


= ( v o +Ax s) 


t 


l\f ( A %) 2 + ( Ax 6 ) ' 


$ 


(v. .4x s ) ; 




Ax c 


neglecting 2 nd order terms and higher. Or, 


A v 0 - AX S 


( 9 ) 


The Error in Flight Path Angle at Insertion 
The error in flight path angle at insertion is 


Ay 0 = y 0 (X + Ax) 


, V v p 

sin 1 ncrro 


or 


(r 0 -‘)ax 2 + (v 0 -‘)ax 4 


ar„ 


(r 0 -)AX 2 + (v 0 -‘)AX, 


( 10 ) 


neglecting 2 nd order terms and higher and noting that the flight path angle is 
zero for a nominally circular orbit. 
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The Error in Energy Per Unit Mass 
The error in energy per unit mass is 

Ac 3 = c 3 (X+Ax) - c 3 (X) 





( 11 ) 


neglecting 2 nd order terms and higher. 
Semi- major Axis Error 


For ease in analysis let k - (r Q v 0 2 //x)* Then, the error in semi- major 
axis may be written 


Aa - a(X + Ax) - a(X) 


( 


r o + Ar o \ ( r o \ 

2 _ ( A + AV ) / \2- k) 


7 



A N +Ar <A f (k A^\ (k A^V 

Aa " \ 2 / y V2 + iy + \t + ~2j + 

“ (t) i 1 + T + fe 


(ft) + (ft) + (( 2 _\) 2 ) ^ “ (ft) + hi i her order t 


erms 


But, for a circular orbit k ~ 1. Hence 


Aa n Ar Q + (r Q ) AA. + higher order terms 


But, 


A \ = (2v Q ^Avq + (r Q Ar 0 + higher order terms 


Thus we can write 


Aa % 2Axj + (2r 0 v 0 _1 ) Ax 5 


or 


AA % 2AXj + ( 2r 0 v 0 ” 1 ) AX S 


( 12 ) 


The Angle Between the Nominal and the Perturbed 
Position Vector at Insertion 


The Angle between R and R p is 


p -- p{\. \) 
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The Error in Eccentricity 

The eccentricity e is a function of r Q , v Q , and y Q 



But for a circular orbit e - 0. Therefore, the error in eccentricity Ae is 

Ae = e(r Q + Ar 0 , v Q + Av 0 , y Q +Ay 0 ) 


9 


or 


AE - e(r 0 +.AR 0 , v 0 +AV 0 , y 0 +Ar 0 ) 


(14) 


The Error in Perigee 

The perigee radius r p is a function of r 0 , v 0 , and 7 0 



For a circular orbit r p = r Q . Therefore, the error in perigee is 

Ar p = r p( r o + A, V v o + Av o • n + A ^o) * r o 


AR P = r p ( r o ,AR o. v„ *AV 0 , y 0 +Ar 0 ) - r 0 
The Error in Apogee 

The apogee radius r a is a function of r Q , v 0 , and y Q 


(15) 
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For a circular orbit r 


r 0 . Therefore, the error in apogee is 


Af a = r a( r 0 + Ar 0- V 0 +AV 0’ >0 + ^o) “ F 0 

or 

^a = r a ( f 0 + AR 0’ V 0 +AV 0’ ^0 + ^o) “ r 0 ( 16 > 


4. THE VARIANCES AND COVARIANCES BETWEEN THE ERRORS 
IN INSERTION RADIUS, SPEED AND FLIGHT PATH ANGLE* 

Using (8), (9), and (10) we may easily obtain expressions for the variances 
of AR 0 , AV 0 , and AT 0 


°A 2 r 0 - ff A 2 ,, (17) 

^ ~ ct A 2 x s (18) 

CT Ar 0 * ( r o“ 2 ) CT A 2 x 2 + (i(T 2 )^ 4 + 2 ( r o v 0 ) _1 ct Ax 2 Ax 4 (19) 


The covariances are 


CTAr o Av 0 


°Ax j Ax j 


( 20 ) 


*If AXj and AXj are random variables, cr^Ax; _ PAx; Axj cr Axj^Axj is the covariance between 
AXj and AXj , PAx, Ax- * s £ ^ e coe Hicient of correlation between them, and cr a x . and cr^ x . 
are the standard deviations of AXj and AXj respectively. 
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CT A V o Ay o 




5 


( 21 ) 


("o' 1 ) CT A 


X 2 Ax 5 


(*.'*) 


a 


Ax Ax 

4 


r 0 Ay 0 - ( r o‘) a Ax,Ax 2 + ( v o _1 ) ^AxjA, 


( 22 ) 


5. THE CUMULATIVE DISTRIBUTION FUNCTIONS OF THE 
ERRORS IN THE PARAMETERS* 

The Cumulative Distribution Functions of the Errors in Insertion Radius , 
Speed, and Flight Path Angle, the Error in Energy Per Unit Mass, and 
the Error in Semi-Major Axis 

From Equations (8), (9), (10), (11) and (12) we see that AR Q , AV 0 , AT 0 , AC 3 , 
and AA (the errors in insertion radius, speed, and flight path angle, the error in 
energy per unit mass, and the error in semi-major axis) are linear combinations 
of the components of AX, and are therefore normally distributed. 

The Cumulative Distribution Func tio ns of the Errors in Eccentricity , 

Perigee Radius, and Apogee Radius 

From Equations (14), (15), and (16) 


AE 

= e ( r o +AR 0’ v o + AV 0’ 

^0 + ^o) 

(14) 


r p ( r o 4 AR 0’ v o + A V 

7o + Ar o) - r o 

(15) 


"a ("0 + ^0- V 0 + AV 0’ 

7o + ^o) - r o 

(16) 


‘The cumulative distribution function F F(At]) of the random variable AH is a function of At] 
and gives the probability that AH is less than or equal to At]. Mathematically, we write 

r At; 

F r F(Atj) : Pr(AH < A 7]) J _ <D f(x) dx 

where f(x) is the probability density function of AH . 
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we see that AE, AR p , and AR a are nonlinear functions of the components of AX, 
and, therefore, are not normally distributed. However, the probability distribu- 
tions of these variables can be calculated. 

Since AR 0 , AV 0 , and AT 0 have a trivariate normal distribution, they can be 
expressed as linear functions of three uncorrelated normal random variables 



^0 S 31 Y 1 + 3 32 Y 2 + a 33 Y 3 

And, by the inverse linear transformation, we can express the Y. 's as func- 
tions of AR 0 , AV 0 , AT 0 

Y 1 “ b 11^0 + b 12^ V 0 + b 13^0 

Yj = b 21 AR 0 * b 22 &v 0 + b 23 ar 0 < 24) 

Y 3 ~ b 31 + b 32^ Y 0 + b 3 3 ^0 

Thus, the variances of the Yj 's 'sj can be calculated as functions of the 
variances and covariances between’ the correlated variables using (24). Since 
uncorrelated normal random variables are independent, the joint density func- 
tion of the Y. 's is given by the product of their marginal frequency functions. 
This is the essence of the technique. 

Each probability density function for Yj is approximated by a discrete dis- 
tribution having probability mass points pj ^ over a range of values y ■ extending 
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from -La to La and spaced, (L/N) a apart.* Thus, 

y i y i y i 



(i = 1, 2, 3; k = 1, 2, •• •, 2N + 1) 


where 



q iu = ( 2 ™r 2 ,) 






for (i = 1, 2, 3; k = 1, 2, 3,---, 2N+ 1) 


and 


S. = 


2N+ 1 

L": 


k= 1 


The larger the integer N for a fixed L, the better the approximation in (25) . 


*a y. is the standard deviation of Yj. 


14 



Since there are ( 2N + 1) possible values for each density function approxi- 
mation, there will be a total of ( 2N + l) 3 possible combinations of values for the 
joint density approximation to the YAs. Further, since the YAs are independent, 
the probability of occurrence of p jkl of each combination of values for y^. 
y 2 k ’ y 3 * . ( j . k . 1 = 1, 2, 3, • • • , 2N+l)is given by the product J 


Pjkl - Pi j P 2k P 3l (26) 

For each set of values for Y 1 , Y 2 , and Y 3 , there is a corresponding set of 
values for AR 0 , AV 0 , and AT 0 by the equations in (23), and a value for one of the 
variables AE, AR p , and AR a (by Equations (14), (15), and (16)) having the proba- 
bility p. kl ( j , k , 1 = 1, 2, 3, 2N+1). Thus, the probability density func- 

tion approximation for AE, AR p , and AR a w iH each have ( 2N + l) 3 possible values. 

For purposes of machine computation, the range of each random variable 
AE, AR p , and AR a can be divided into a number of mutually exclusive intervals. 

For all values of the variable falling within one of these intervals, the associated 
probability as given by Equation (26) can be summed since all of the ( 2N + l) 3 
possible combinations of values for AR 0 , AV 0 , and AT 0 are mutually exclusive. 

In this way we obtain an approximation for the probability that the random vari- 
able (AE, AR p , AR a ) assumes values within the particular interval, and thus an 
approximation for the probability density function. 

The cumulative distribution function is obtained by summing the probability 
density function. 

The Cumulative Distribution Function of B 


By Equation (13) 


B a r 0 -> j/(AX 2 )2 + (AX,) 2 


we see that B is not a linear function of the components of AX and therefore, not 
normally distributed. The cumulative distribution function of Bean be obtained 
in similar fashion as for AE, AR p , and AR a . However, in this case B is a function 
of only two components of AX, AX 2 and AX 3 . 
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6. AN EXAMPLE WITH APPLICATION TO A NOMINAL 
APOLLO EARTH PARKING ORBIT 

In order to illustrate the statistical technique described in Section 5, we will 
consider a nominal Apollo earth parking orbit— a nominally circular orbit of 
100 n.mi. 

Reference 3 describes the dispersions in position and velocity of the Saturn 
V Launch Vehicle at insertion due to 30 navigation parameters. Reference 4 
describes how the covariance matrix of insertion errors shown in Table 1 and 
obtained from Reference 5 was constructed from the uncertainties (standard 
deviations) of these navigation parameters. The matrix, therefore, expresses 
the deviation of the onboard estimate from the nominal insertion condition due 
to errors in the onboard navigation during the launch phase. It can be considered 
typical of the performance of the Saturn V Launch Vehicle for a nominal Apollo 
Mission at insertion into a near-earth nominally circular parking orbit of 100 
n.mi. as of the date of this analysis (Reference 6). The matrix is given in the 
U° V° W° Coordinate System (Figure 1). The elements along the main diagonal 
are the variances of AX ; (the components of AX) expressed in units of (ft) 2 and 
(ft/sec) 2 . The off-diagonal elements are the covariances between AX. and 
AX j (i / j) and are expressed in units of (ft) 2 , (ft/sec) 2 , and (ft) 2 /sec. The 
normalized covariance matrix or correlation matrix corresponding to the matrix 
in Table 1 is shown in Table 2. It is also given in the U° V° W° Coordinate System. 
However, the elements along the main diagonal are unit variances while the off- 
diagonal elements are the coefficients of correlation between AX. and AX^ ( i f j ) . 

It is necessary to mention at this point that we are assuming throughout this 
analysis that AX has a multivariate normal distribution with zero mean and the 
covariance matrix shown in Table 1. 


Table 1 

Expected Saturn V Insertion Covariance Matrix* ( U° V° W° Coordinate System) 


977736.00 -745996.00 


Symmetric 


743820.00 


-3162.1120 

1493.9920 

1180016.0 


3935.6520 

-3123.1360 

-18.365520 

16.197720 


-2060.5760 

1961.3040 

5.3787600 

-8.6568800 

5.2908400 


-15.106120 

8.9479200 

3499.8720 

-0.0769048 

0.02838016 

11.023800 


2 

•The elements along the main diagonal are the variances of AX i expressed in units of (ft) and 
(ft/sec) 2 , the off-diagonal elements are the covariances between the AX*, expressed in units of 
(ft) 2 , (ft/sec) 2 , and (ft) 2 /sec. 
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One may ask, "What actually does a covariance matrix of insertion errors 
represent?" And, "What is the usefulness of such a matrix?" 

The covariance matrix is actually a measure of the uncertainty in the state 
vector (position and velocity) at insertion due to uncertainties in the actual values 
of navigation parameters. If the guidance and control system is working properly, 
we could still expect this much variation or dispersion in the position and velocity 
of the vehicle at insertion. 

The usefulness of such a matrix can best be seen in its providing us with a 
tool for determining expected variations and bounds (along with their associated 
probabilities under the Gaussian assumption) of some of the parameters of the 
orbit. 

Using the expressions for the variances and covariances between AR 0 , AV Q , 
and r o in Equations (17) through (22) as well as the covariance matrix in Table 1, 
the covariance matrix of AR Q , AV 0 , and AT 0 can be obtained, and is shown in 
Table 3. Table 4 shows a normalized matrix where the elements along the main 
diagonal are the standard deviations of AR AV 0 , and AT 0> and the off-diagonal 
elements are the coefficients of correlation between these random variables. 

Table 3 

Covariance Matrix S 0 ofAR 0 , AV 0 , and AT 0 * 

0.02644932 -0.33891053 -0.00100201 

5.29084000 0.01588871 

Symmetric (0.49300523)10" 4 

*The elements along the main diagonal are the vari- 
ances of ARo, AVq, and APg, expressed respec- 
tively in units of (n.mi.) 2 , (ft/sec) 2 , and (deg) 2 . 

The off-diagonal elements are the covariances be- 
tween ARq, AVq, and AFq, expressed in units of 
(n.mi.)(ft/sec), (n.mi.)(deg), and (ft)(deg)(sec). 


The cumulative distribution functions for some of the errors in the parking 
orbit are shown in Figures 2 through 10, and were obtained as described in 
Section 5. 

From these figures we can see that the angular error in the position vector 
at insertion (Figure 7), the error in eccentricity (Figure 8), the perigee error 
(Figure 9), and the apogee error (Figure 10), are non-Gaussian and do not have 
zero means. 


18 



F, (Ar 0 ) 


Table 4 


Normalized Covariance Matrix 
of AR Q , AV 0 , Ar 0 * 

0.16263246 -0.90597355 -0.87748401 

2.30018258 0.98378570 

Symmetric (0.70214331)10" 

•The elements along the main diagonal are the 
standard deviations of ARq, AVq, and AFo, ex- 
pressed in (n.mi.), (ft/sec), and (deg). The off- 
diagonal elements are the coefficients of correlation 
between ARq, AVo, and AFq, 



(n.mi .) 


Figure 2. The Cumulative Distribution Function of the Error in Insertion Radius ARq 
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F igure 5. 


The Cumulative Distribution Function of the Error in Energy Per Unit Mass AC3 
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Figure 6. The Cumulative Distribution Function of the Error in Semi-major Axis AA 
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Figure 7. The Cumulative Distribution Function of the Angle Between the Actual Position Vector 

and the Nominal Position Vector at Insertion B 
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Figure 8. 


The Cumulative Distribution Function of the 


Error in Eccentricy 
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Figure 9. The Cumulative Distribution Function of the Error in Perigee Radius AR p 
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Figure 10. The Cumulative Distribution Function of the Error in Apogee Radius AR a 
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Table 5 


Lower and Upper 99 .5% Probability Points for the Errors 
in the Parameters of the Orbit 


Error 

Lower 99.5% Point 

Upper 99.5% Point 


-0.7 km (-0.4 n.mi.) 

+0.7 km (+0.4 n.mi.) 

AV o 

-1.8 m/sec (-6.5 ft/sec) 

+1.8 m/sec (+6.5 ft/sec) 


-0.018° (-0.3 mrad) 

+0.018° (+0.3 mrad) 

ac 3 

-0.016 (km/sec) 2 

+0.016 (km/sec) 2 

AA 

-1.7 km (-0.9 n.mi.) 

+ 1.7 km (+0.9 n.mi.) 

B 

+0.0002° 

+0.0085° 

AE 

+0.0001 (dimensionless) 

+0.00048 (dimensionless) 


-4.7 km (-2.5 n.mi.) 

+0.2 km (+0.1 n.mi.) 

A *a 

-0.2 km (-0.1 n.mi.) 

+4.5 km (2.4 n.mi.) 


Table 5 summarizes the results of this analysis in tabular form, showing 
both lower and upper 99.5% probability points for the errors in the parameters. 
Thus, there is a 99% probability that the errors in these parameters will lie 
between these indicated values. 

It is significant to note that there is a 99.5% probability that perigee error 
will exceed -2.5 n.mi. or that the actual perigee height will exceed 97.5 n.mi. 
for the nominal 100 n.mi. circular parking orbit. Likewise, there is a 99.5% 
probability that the actual apogee height will be less than 102.4 n.mi. 

From the above we see that the actual parking orbit should be very close to 
the nominal circular orbit of 100 n.mi. This is also evident from inspection of 
Figure 8 where we see there is a 99.5% probability that the error in eccentricity, 
and hence the eccentricity, will be less than 0.00048. 

Figure 11 shows a comparison between the cumulative distribution function 
of perigee error using Equation (15) and a cumulative normal distribution func- 
tion having the same mean and standard deviation. The difference between the 
curves is significant. For example, using the normal distribution, there is a 
99.5% probability that the perigee error will exceed -1.7 n.mi. and a 99.5% 
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PROBABILITY THAT AR p *Ar p 
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Figure 11. Comparison of the Cumulative Distribution Function of the Error in Perigee with a 
Cumulative Normal Distribution Having the Same Mean and Standard Deviation 


probability that the error will be less than 1.1 n.mi. From the calculated curve 
the corresponding points are -2.5 n.mi and 0.1 n.mi. 

Combining Tracking Errors and Insertion Errors 

As a further example, we will show how insertion ship tracking errors (cal- 
culated minus actual) can be combined with insertion errors (actual minus 
nominal) to answer such questions as, "For a nominally circular parking orbit 
of 100 n.mi., what is the probability that the calculated perigee height will exceed 
a given amount? " 

At insertion of the Apollo spacecraft, the following parameters: 

radius r t 

speed Vj 

flight path angle 7 1 
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will be determined from one minute of ship’s tracking data. In Reference 1 it is 
required that these parameters be determined to the following accuracies: 

3cr = 2.4 n.mi. (4.44 km) 

r i 

3 a v = 16 ft/sec (4.87 m/sec) 

3a = 0.16° (2.79 mrad) 



be an error vector representing the error or deviation of the actual insertion 
radius, speed, and flight path angle from the nominal , and let 




ARj 

AV X 


(28) 


be an error vector representing the deviation of the measured or calculated 
insertion radius, speed, and flight path angle from the actual . 

Then 


AZ - AZ Q + AZj (29) 

is an error vector representing the deviation of the measured or calculated in- 
sertion radius, speed, and flight path angle from the nominal . 

For this example, we will assume that AZ 0 and AZj are independent error 
vectors each with a trivariate normal distribution, having zero mean vectors and 
covariance matrices Xq and X 1 respectively. Then AZhas a trivariate normal 
distribution with a zero mean vector and covariance matrix X 

X X Q + Xj (30) 
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Since it is practically impossible to determine the actual coefficients of 
correlation between the errors AR X , AV X , and AT t , the analysis has been carried 
out for coefficients of correlation of 0.0 (uncorrelated) and ±0.9. That is, all 
possible combinations of signs have been considered with the expection of those 
for which the covariance matrix 2 j becomes singular. The possible combinations 
are: 


^Ar, Av, 
+ 0.9 

+ 0.9 
-0.9 
- 0.9 


l °Ar 1 A-y, 
+ 0.9 

- 0.9 
+ 0.9 

- 0.9 


^Avj Ay, 

+ 0.9 

- 0.9 

- 0.9 
+ 0.9 


(31) 


The matrix 2j is shown in Table 6 for the case p^ T 
^Av,Ay, = 0,9 ‘ 


/0 Ar [ Ay, 


Table 6 

Covariance Matrix 2j of 

ARj , AVj , andAfj* 


0.64 

3.84 

0.038400 


28.4444 

0.25600 

Symmetric 


0.002844 


*The elements down the main diagonal are the 
variances of ARj, AV i, and AFi expressed in 
units of (n.mi.)2 ; (ft/sec)2, and (deg)2. The off- 
diagonal elements are the covariances between 
ARj, AV i, and Ar j, expressed in units of 
(n.mi.)(ft)/sec, (n.mi.)(deg), and (ft)(deg)/sec. 


The matrix 2 obtained by adding 2 0 (Table 3) and 2 j (Table 6) is shown in 
Table 7. Similar matrices can be obtained using other coefficients of correla- 
tion indicated in (31) and A = PA fj = PAv, Ay, = 0 0 - 

Using these matrices and the method outlined in Section 5, the probability 
distribution of the perigee error can be calculated. 
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PROBABILITY THAT CALCULATED PERIGEE EXCEEDS r, 


Table 7 


Covariance Matrix X of (AR 0 + AR^ , 
(AV 0 +AVJ , and (AT 0 +A rj* 

0.66644932 3.50108947 0.03739799 

33.73524000 0.27188871 

Symmetric 0.00289330 

•The elements along the main diagonal are the 
variances of (ARq + AR j), (AVo + AV i), and 
(AP o + AP i), respectively, expressed in units of 
(n.mi.) 2 , (ft/sec) 2 , and (deg) 2 . The off-diagonal 
elements are the covariances between these ran- 
dom variables, expressed in units of (n.mi.)(ft)/ 
sec, (n.mi.) deg, and (ft)(deg)/sec. 



(n.mi . ) 
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Figure 12. Probability That Calculated Perigee Height Exceeds r p ^ 
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Figure 12 shows the cumulative distribution function of the calculated perigee 
height for a 100 n.mi. parking orbit and several combinations of coefficients of 
correlation between ARj , £V 1 , and ATj . 


In order to compare these curves, we note that if p^ r = + 0.9, and 
p^r = pp v = ±0.9, there is a 90% probability that the calculated perigee 
height is greater than 90.8 n.mi. If the errors are uncorrelated, there is a 90% 
probability that the calculated perigee will exceed 92.3 n.mi. If p Ar p v = - 0.9, 
p^ r = ± 0.9, and p Av = + 0.9, there is a 90% probability that the calcu- 
lated perigee height will exceed 93.7 n.mi. 


CONCLUSIONS 


Distributions 


A statistical technique has been described for analyzing the effects of in- 
sertion errors (deviations of actual values from nominal values) on the param- 
eters of near-earth circular orbits. The technique assumes that the error in 
the state vector (position and velocity) of the vehicle at insertion has a multi- 
variate normal distribution with zero mean and a given covariance matrix. 
Probability distributions for some of the errors in the parameters are obtained 
numerically using Keplerian two-body equations of motion. Results of the 
analysis indicate: 

(1) Perigee error, apogee error, the error in eccentricity, and the angular 
error in the vehicle's position vector at insertion are non-Gaussian 
with non- zero means. 

The Apollo Parking Orbit 

The technique has been applied to a nominal Apollo earth parking orbit using 
an expected covariance matrix typical of the performance of the Saturn V Launch 
Vehicle at insertion into a near-earth nominally circular parking orbit of 100 
n.mi. Results of the analysis indicate that: 

(2) The actual near-earth parking orbit will be close to the nominally cir- 
cular 100 n.mi. orbit. For example, there is a 99.5% probability that 
the actual perigee height will exceed 97.5 n.mi., and a 99.5% probability 
that the actual apogee height will be less than 102.4 n.mi. 
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Combining Tracking Errors and Insertion Errors 


As a further example, insertion ship tracking errors and insertion errors 
have been combined. The analysis indicates that for 3 sigma values of 


3 cta = 2.4 n.mi. (4.44 km) 

nr j 

3cr^ v = 16 ft/sec (4.87 m/sec) 
3 ct A 7i = °- 16 ° (2.79 mrad ) 


for the insertion ship’s tracking errors (calculated minus actual) in insertion 
radius r j , speed Vj, and flight path angle y l , and for coefficients of correlation 
of +0.9 between these errors: 

(3) There is a 90% probability that the calculated perigee height will exceed 
91 n.mi. for the nominal 100 n.mi. parking orbit. 
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